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Abstract 


Logical RS flip-flop circuits are investigated once again in the context of discrete 
planar dynamical systems, but this time starting with simple bilinear (minimal) com¬ 
ponent models based on fundamental principles. The dynamics of the minimal model 
is described in detail, and shown to exhibit some of the expected properties, but not 
the chaotic regimes typically found in simulations of physical realizations of chaotic RS 
flip-flop circuits. Any physical realization of a chaotic logical circuit must necessarily 
involve small perturbations - usually with quite large or even nonexisting derivatives 
- and possibly some symmetry-breaking. Therefore, perturbed forms of the minimal 
model are also analyzed in considerable detail. It is proved that perturbed minimal 
models can exhibit chaotic regimes, sometimes associated with chaotic strange attrac¬ 
tors, as well as some of the bifurcation features present in several more elaborate and 
less fundamentally grounded dynamical models that have been investigated in the re¬ 
cent literature. Validation of the approach developed is provided by some comparisons 
with (mainly simulated) dynamical results obtained from more traditional investiga¬ 
tions. 
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Since the advent of electronic computing devices, flip-flops have played a 
major role in processing arithmetic operations. The first flip-flop to be designed 
was the RS flip-flop (RSFF). R has the ability to store bits, and hence was used 
in adders and primitive electronic calculators. The RSFF was even used in the 
video game “pong”. The only logical flaw of this circuit is the ambiguity of the 
output when the inputs are both high voltage. In order to rectify this flaw other 
flip-flops were designed. In the 1980s Leon Chua designed the simplest circuit 
to exhibit chaos. In more recent years, Chua’s circuit and the flaw of the RSFF 
have been exploited to design chaotic flip-flops. These chaotic flip-flops have 
the potential to be employed in random number generation, encryption, and 
fault tolerance. However, in order to exploit their properties they need to be 
studied further. Since experiments with large systems become difficult, tractable 
mathematical models that are amenable to analysis via the tools of dynamical 
systems theory are of particular value. The model needs to be simple enough to 
simulate with ease and agree qualitatively with experiments for small systems. 
We model the RSFF circuit as a discrete dynamical system - called the minimal 
model - focusing on the ambiguous output case. The minimal model is perturbed 
in various ways to reproduce qualitative behavior observed in experiments for 
the chaotic RSFF, and the minimal approach can be readily generalized to handle 
other logical circuits. 


1 Introduction 


Logical circuits are constructed using logic gates representing propositional connectives such 
as AND (conjunction), OR (disjunction) and their negation. They, or more precisely physical 
approximations of ideal logical circuits, have important applications in multiplexers, regis¬ 
ters, pseudo-random number generators, quantum network modeling and, in fact, virtually 
every microprocessor (see, e.g. 0). Consequently, having efficient methods for analyzing and 
predicting their behavior, such as by continuous or discrete dynamical systems models is of 
great value in design, analysis and evaluation. 

A well-known example with many applications {e.g. in quantum networks [31]) is the RS 
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Flip-Flop circuit in Fig.l, which is a feedback circuit comprising two NOR gates (which are 
negations of OR gates). Note that an OR gate has an output of 1 when at least one of the 
inputs is 1 and an output of 0 when both inputs are zero, with 1 denoting true and 0 false. 
An ideal RS flip-flop circuit [RSFF circuit) is a logical feedback circuit represented in Fig. 
1, with input/output behavior described in Table 1, which shows the set (S) and reset (R) 
inputs for the circuit consisting of two NOR gates with outputs Q and Q'. The input to the 
output, denoted by {Qn, Q'fl) {Qn+i, Q'n+i) be regarded as the action of a map from 
the plane := {{x,y) : x, y G M} into itself, where M denotes the real numbers and in the 
ideal case, the coordinates assume the binary values {0,1}. This suggests that the behavior of 
the RSFF circuit might possibly be effectively modeled by the iterates of a planar map, which 
comprise a discrete dynamical system (DDS), and its perturbations derived directly - albeit 
in a simplihed manner - from a faithful interpretation of the circuit itself. And this is our 
focus in the sequel, which should be contrasted with our prior investigation j^. The model 
to be studied in what follows, fundamentally grounded as it is in the basic principles of ideal 
logic circuit theory, is substantially more compelling and useful than the one investigated 
in [2], which was primarily an ad hoc construct designed to mimic the known behavior of 
RSFF realizations. 



The binary input/output behavior, with 0 and 1 representing false and true, respectively, is 
given in the following table. 


A 

R 

Si:=Q 

Ri := P 

1 

0 

1 

0 

0 

1 

0 

1 

1 

1 

0 

0 

0 

0 

1 or 0 

0 or 1 
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Table 1. Binary input/output of R-S flip-flop circuit 


From the DDS perspective, our first goal is to construct the simplest map of the plane 
(based on bilinear representations of the NOR gates) that models the logical properties of 
the RSFF circuit, with iterates that exhibit most of the interesting properties that follow 
from basic analysis or have been observed in the dynamics of physical realizations and their 
mathematical models. This begs the question of how the dynamics of the planar map models 
are to be compared with that obtained directly from measurements of physical realizations 
of the RSFF circuit and its flip-flop relatives and the analysis and simulations of the usual 
mathematical models, which we shall endeavor to address in what follows. 

The usual form of the mathematical models is traceable back to the pioneering work 
of Moser [21] ; namely, associated three-dimensional systems of piecewise-smooth, hrst- 
order, nonlinear autonomous ordinary differential equations (ODEs) obtained from applying 
Kirchhoff’s laws to the realizations. These realizations typically comprise such elements as 
capacitors, inductors and nonlinear resistors and exhibit highly oscillatory, very unstable and 
even chaotic dynamics {metastable operation), as experimentally observed in such studies as 
pEl [I7112I], where tunnel diodes of the type used in Chua’s circuit (see aElElES] ) 
are the key ingredients in the construction of the nonlinear resistors. There are several 
connections between the solutions of the model ODEs and iterates of maps that can be used 
for dynamical comparisons, among which are the following: as observed by Hamill et al. 
[I3] . the autonomous nature of the logical circuit equations allows dynamical analysis via 
the iterates (snapshots) of a hxed time map, which can be reduced to a planar map in special 
cases as shown in Kacprzak and Albicki [16] and Kacprzak mi. In addition, the application 
of a standard explicit one-step integration method is tantamount to the iteration of the map 
describing the scheme, thus enabling the (approximate) reduction from a continuous to a 
discrete dynamical system as shown e.g. in Danca [8]. Of course, there is the well-known 
method of employing Poincare sections to analyze three-dimensional continuous dynamical 
systems using two-dimensional discrete dynamical systems, which has been employed in 
numerous investigations of logical circuit realizations such as Cafagna and Grassi [1] , Murali 
et al. [26], Okazaki et al. mi and Ruzbehani et al. |31j . 

As it turns out, logical circuits can be realized using the famous circuit of Chua and its 
generalizations [a [7], which have depended heavily powerful tools such as Poincare maps, 
Melnikov functions and normal forms from the modern theory of dynamical systems and 
bifurcation theory (see 0 El tta EDI ESI El ES]) for their analysis. There is also another 
interesting connection between realization of logical circuits and nonlinear maps that might 
afford an opportunity for comparisons with the results obtained from our two-dimensional 
DDS models; namely, the approach in the work of Ditto et al. m, which features the notion 
of reconhgurable logic gates comprising connected NOR and NAND gates constructed using 


4 


1-dimensional DDS and associated thresholds. 


Our investigation begins in Section II, where we define a (minimal) planar map model of 
an ideal RSFF circuit - derived directly from the logical circuit - that plays a foundational 
role throughout the sequel. Moreover, we derive some basic properties of the minimal map 
concerning such things as smoothness and the existence and description of an inverse. This 
is followed in Section III with a more thorough analysis of the fixed points of the minimal 
map - including a local stability analysis and an analysis of stable and unstable manifolds. 
As a result of this more detailed investigation, we find that the dynamics is quite regular 
when the domain of the map is appropriately restricted. Next, in Section IV, after observing 
that the minimal map is C^-structurally stable, we prove that (1-dimensional) chaos can 
be generated by arbitrarily small perturbations of the map. In particular, it is shown 
that a tent map can be embedded in a (1-dimensional) stable manifold of a fixed point by 
such a perturbation, thereby inducing chaotic tent dynamics. We then prove the existence 
of several types of more substantial (2-dimensional) chaotic regimes for arbitrarily small 
perturbations in Section V. In particular, we show that arbitrarily small perturbations 
can be fashioned to produce transverse intersections in homoclinic orbits and heteroclinic 
cycles that generate chaos, horseshoe chaos, multihorseshoe strange chaotic attractors, and 
Neimark-Sacker bifurcations. In Section VI, we illustrate our theorems by making liberal 
use of numerical simulations of our perturbed models. Finally, in Section VII, we briefly 
summarize our results, describe some interesting areas of application and discuss plans for 
related future research. 


2 A Minimal Discrete Dynamical System Model 


To obtain the minimal DDS model for the RSFF circuit, we start with the simplest continuous 
extension of the input/output map of the NOR gate to the unit square in the plane, which 
is clearly given by the function M P (c M^) —)■ / (c M) dehned as 

M{x,y) ■.= {I-x){l-y), (1) 

where P := I x I := [0,1] x [0,1] C is the unit square and x and y are the inputs. 
Observe that this quadratic function is bilinear in the variables (1 — x) and (1 — y), and 
behaves precisely like the pure logical NOR gate when x,y & {0,1}. Using the functional 
representation ([^ for each of the (NOR gate) components of the RSFF circuit, we obtain 
the continuous extension for the RSFF, 


^( 1 -^) w . ^(1 - g ) 

R + S -RS' ^ ' R + S-RS' 


( 2 ) 
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Since the threshold outputs are fed back in a similar manner to the RSFF ouputs, we 
make the naive assumption that the logical behavior is preserved and we readily obtain its 
minimal extension to the unit square in the form of a map F : P ^ P dehned as 


F = {^,ri) : P (c M^) ^ P (c M^) , 

where the coordinate functions are 




?/(l - x) 

X + y — xy' 


7]{x,y) 


xjl - y) 

X + y — xy^ 


(3) 

(4) 


with the X- ?/-coordinates playing the roles of the previous threshold ouput, and ^ - rj coor¬ 
dinates playing the roles of the current threshold ouput. 


This map generates a discrete (semi-) dynamical system in terms of its forward iterates 
determined by n-fold compositions of the map with itself, denoted as F”, where n G N, the 
set of natural numbers. We shall employ the usual notation and dehnitions for this system; 
for example, the positive semiorbit of a point p G which we denote as 0+{p), is simply 
dehned as 0+(p) := {F^{p) : n G Z, n > 0}, and all other relevant dehnitions are standard 
(c/. [121 dSl HSl 132] )• The minimal map (|^ is real-analytic (C"^) on P except at the origin, 
where it is not even well-dehned. Consequently, we should really consider F to be dehned 
on P \ {(0, 0)} and actually on X := P \ {(0, 0), (1,1)} if we wish to avoid the origin for 
all forward iterates of F. 


2.1 Basic analytical properties of the minimal model 

First we take note of some of the analytical properties of the map ([^, which are listed in 
what follows. The proofs of all these results are straightforward and left to the reader. 

(Al) The map F G {X). 

(A2) F is Z 2 -symmetric in the sense that 

R o F = F o R 

for the rehection R in the line y = x. 

(A3) F (X) = F := {(e, p) G M2 ; 0 < e, p and e + P < 1} U {(1, 0), (0,1)}. 
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(A4) F ({(x, 0) : 0 < a; < 1}) = (0,1), F ({(0, y) : 0 < y < 1}) = (1, 0), and if p(a) := {(x, ax) : 
0 < x, a} n X is the ray through the origin sans the origin in X, then F (p(a)) C 
{(■C) ~ 1) + : 0 < .^} n F. Moreover, as {x, y) —)■ (0, 0) along the ray p(a), F(a;, y) 

converges to the point of intersection of the line dehned hy x + y = 1 with the ray 
p (1/a). This shows just how singular the formulas (3) are at the origin. 

(A5) The derivative (matrix) for F on X is 

F'{x,y) = {x + y-xy)~^ (5) 

with determinant 

det F' {x, y) = xy {x + y — xy)~^ , (6) 

which, not surprisingly in view of (A3), shows that the implicit function theorem cannot 
guarantee the existence of a local smooth inverse along the x- and p-axes. 


(A6) In fact, the inverse of F, where it exists, is given as 




V 1-^ ’ i-e A 


which is clearly in (Y \ {(1, 0), (0,1)}). 


(7) 


3 Dynamics of the Minimal Model 


We shall analyze the deeper dynamical aspects of the model map (2)-(3) for various parameter 
ranges in the sequel, but hrst we dispose of some of the more elementary properties that 
follow directly from the dehnition and (Al)-(A6), leaving the simple proofs once again to 
the reader. 


(Dl) If we restrict F to X := {{x,y) : 0 < x,y and a; + p < 1}, and denote this restriction 
by F, it determines a full dynamical system dehned as 

|A:nez|, 

which, for example, allows the dehnition of the full orbit of a point p G F as 

0{p) ■= I A(p) : n e z|. 

(D2) The line y = x is F-invariant, while the x- and p-axes are F^-invariant. 

(D3) Both of the points (1, 0) and (0,1) are hxed points of F^. 
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3.1 Analysis of the fixed and periodic points 


The properties of the hxed and periodic points of onr model map shall be delineated in a series 
of lemmas, which follow directly from the resnlts in the preceding sections and fnndamental 
dynamical systems theory (as in [T^[T8| l28 | 133]). Our first result is the following, which has 
a simple proof that we leave to the reader. 


Lemma 1. The only fixed point of F in X is 
p* = = I — 

which is a saddle point with eigenvalues 

'3 - Vs' 


(1,1) = 0.38197(1,1), 


A. = 




'1 + ^5' 


This fixed point has linear stable and stable manifolds given as 

^tin (P*) = {(a:, x) : X e M} and IT® (p*) = {(x, x) : 0 < x < 1}, 
and the linear unstable manifold 

ip*) = -X + 2xd : X e M}. 


Next, we analyze the unstable manifold of p* in some detail. For this purpose, it is 
convenient to introduce a change of variables linked to the symmetry of F expressed in (A2); 
namely. 


T{x,y) = {u,v) 


/x + p - 2x* 

V 71 


-X + y\ 

72 )' 


with inverse 


T \u,v) = (x,p) 



u — V 



+ 


u + v\ 

~w)' 


Clearly, T is a translation of the origin to the hxed point p* followed by a counterclockwise 
rotation of vr/d. The dehning map for the dynamics in the new coordinates can be readily 
computed to be 


F{u,v) :=T o F oT ^(n, n) = (^V2/R{u,v)^ (^72(1 + 2xl)u + (1 — — u^), —\/2nj , 

( 8 ) 

where 

R{u, v) := 2x*(2 — x*) + 2a/2(1 — x*)m + 

The properties of the unstable manifold of the hxed point may now be described as in the 
following result, which can be proved directly from Lemma 3.1 and (7). 









Lemma 2. The unstable manifold of the fixed point p*, which corresponds to 0 = (0,0) in 
the new uv-coordinates, has the form 

(p,) = (0) = :\v\< 1/v^}, 

where (p is a smooth {C°°) function satisfying the following properties: 


(i) </9(0) = ip'iff) = 0 and ip(y) f as v f l/\/2. 

(a) ip is an even function 
(Hi) ip satisfies the functional equation 

(p{v) = [(p{v) -V ) -i -- (p 


V2 

where k := 1/ [ a /2 (1 + 2x1)] 

S{p:{v),v) := 2 x *(2 — a :*) + 2 a / 2(1 — x^)p){v) + 


(9) 


We note here that (8) can be used to obtain Picard iterate (local) approximations of the 
unique solution via the recursive formula (c/. |14j ) 


( \ //I 2 \ ( ( \2 2A , ‘S'((p„(n),n) 

P>n+l(v) = K I (1 - a:J [p)n{v) -V )+ -- p>r 

A good starting point for these iterates is 

(pi(n) := a/2 ^a/s - 2^ 


-2v 


S{p>n{v),v) 


( 10 ) 


which satishes the hrst two properties of Lemma 3.2, and turns out to be a fairly good 
approximation for the unstable manifold. Using (9), we obtain an even better approximation 
in the form 


(P2{u) = K 


where 


(l-xl) (2(75-2)7 


+ 


a/ 5-2) S'((pi(n),n) 


-2v 


S{ipi{v),v 



( 11 ) 


S ((pi(n),n) = 2 ja:* (1 - X*) + (2(1 - x*)(\/5 - 2) + l) - (Ts - 2) 


\2 2 


which is illustrated in Fig. As a matter of fact, it can be proved that these iterates actually 
converge locally to the smooth solution of (8), but the details, which follow the argument 
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in HI, although straightforward, are a bit too involved to include here. We remark here 
that it is not difficult to prove that a global solution of the unstable manifold equation can 
obtained as follows: 

(p,) = lim„^oo (A'), 

where A' := {(x, 1 — x) : 0 < x < 1. 

Let us now investigate periodic orbits of the dynamical system (2). It is easy to see from 
the definition that ( := {(1, 0), (0,1)} is a 2-cycle in which each point has (least) period 
two. As for any other cyclic behavior, we have the following comprehensive result that 
follows directly from the definition of the dynamical system and properties (Al)-(A6) and 
(D1)-(D3). 

Lemma 3. The 2-cycle C. ■= {(1,0), (0,1)} is the only cycle of F; it is superstable and has 
basin of attraction 


sB(() ;= /2 \ : X e M}. 


Ap proximati on of t h e Unotabl e Meuiif ol d in u,v< o ordinal es 


Ap proMmati on of t ha tMatable Manif ol d in x,y< o ordinal as 




Figure 2: Picard iteration approximation [ip^) in m, x-coordinates and x, ^-coordinates re¬ 
spectively. 


3.2 Summary of the dynamics 


We see from our analysis that the dynamical system (generated by) F defined by (2)-(3) has 
highly oscillatory, but quite regular behavior. There is one fixed point (on the line y = x), 
which attracts everything on the diagonal in X. In addition, there are only two other special 
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points; namely, (1,0) and (0,1), both of which are superstable periodic points of period two 
that comprise the 2-cycle ( := {(1,0), (0,1)}. Moreover, ( attracts everything in X except 
the points along the diagonal. Thus we see that our minimal model nicely replicates the 
highly oscillatory “race” behavior found in physical realizations of RSFF circuits, but none 
of the chaotic dynamics. In fact, the nonwandering set of the minimal dynamical system, 
which actually coincides with the periodic set, has the very simple form 

n (F) = Per (F) = {p4 U C, (12) 

and F is C'^-structurally stable so that the dynamics maintains its regularity for all pertur¬ 
bations that, along with their derivatives, are sufficiently small. A proof of the structural 
stability for our map, which does not quite satisfy the usual hypotheses, can be fashioned 
from a straightforward modihcation of the methods employed in P ESI ESI ED]- Thus, it 
would seem that perturbations capable of generating chaos can be very small, but their 
derivatives need to be quite large or even fail to exist. We shall verify this in the sequel. 


4 Perturbed Minimal Model with One-dimensional Chaos 


In this section we shall show how small perturbations that do not break the reflectional 
symmetry of the minimal model can produce chaos along the diagonal. More precisely, we 
shall give a constructive proof of the following result. 

Theorem 4. There exist eontinuous arbitrarily small perturbations of the minimal model 
map F for the RSFF cireuit that are symmetric with respect to reflections in the invariant 
line y = X and exhibit one-dimensional chaos in their restrictions to the diagonal. 


Proof. The idea of our argument is to embed an arbitrarily small perturbation into F that 
implants known chaotic dynamics along the dynamics without breaking the symmetry. This 
chaotic insert is dehned along the diagonal for any a > 0 using the continuous piecewise- 
linear function 


'lfa{x) 


—X, —2a < X <0 

—2{a—\x — a\), 0 < X < 2a ’ 


(13) 


which is illustrated in Fig. 8a Observe that the composition To- := ifl : [0, 2a] —)■ [0, 2ct] is 


just the 2cT-scaled tent map given by 


rflx) = 2{a -\x- a\), 


which is known to have chaotic dynamics, including periodic orbits of all periods, a dense 
orbit and a Lyapunov exponent of log2 for almost all initial points (see, e.g. [IB])- Now we 
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note from the definitions (2) and (3), (A2) and Lemma 3.1 that the restriction / of F to the 
diagonal is 

m = ( 14 ) 

which has the nniqne (stable) hxed point a:* = (1/2) (3 — \/5), and this is a global attractor 
on the unit interval [0,1]. It is easy to see that for every e > 0, we can choose 0 < Sa < 
min{l/6,e} such that there is a continuous map fa : [0,1] —t [0,1/2] with the following 
properties: (i) fa{x) = tfa^x — x^) for all x G := [x* — 2(j, x* + 2a]; (ii) /o-(x) = /(x) for 
all X G [0,1] \ [x* — 8(7, X* + 4(j]; (hi) fa is strictly decreasing on [0,x* + a] U [x* + 2a, 1]; 
and (iv) |/o-(x) — /(x)| < 4cr for all 0 < x < 1. We observe that Ja is a global attractor for 
fa, but not strange since it has dimension equal one. The one-dimensional chaotic implant 
can be extended to the whole two-dimensional domain of F. More precisely, it is easy to 
see (by simply linearly joining the perturbed diagonal map to F) that there exists for any 
0 < 8cr < min{l/6, e} a continuous map Fa : X ^ X satisfying the following properties: (a) 
the diagonal is Fa- invariant (b) Fa restricted to the diagonal is fa] (c) Fa is symmetric with 
respect to the diagonal; (d) the dynamics of Fa of the diagonal are qualitatively the same as 
that of F; and (e) \Fa{x,y) — F[x,y)\ < 8(j for all {x,y) G X. As a can be made arbitrarily 
small, the proof is complete. □ 


Observe that the perturbation Fa constructed in the above proof is merely continuous and 
piecewise smooth inside a small ball centered at the hxed point x* and smooth (i.e.O^) 
outside of the ball. The perturbation and its chaotic dynamics are illustrated in Fig. 8b 


By smoothing the corners of the construction, the perturbation can be made smooth, which 
leads directly to the following result. 


Corollary 5. The perturbation in Theorem can be chosen so that it is a function 
having the same qualitative dynamics as the function constructed above. 


It is useful to note that the perturbation Fa producing one-dimensional chaos along the 
invariant diagonal is arbitrarily close to the original map F. However, it cannot be 
made arbitrarily enclose to F by virtue of the readily verihed fact that the restriction to 
the diagonal is C^-structurally stable, so we cannot hnd arbitrarily small perturbations in 
the sense that posses chaotic dynamics. As we shall see in the next section, there exist 
arbitrarily small perturbations of F that are C°^ and exhibit (two-dimensional) chaotic 
dynamical regimes. 
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5 Two-dimensional Chaos Induced by Perturbation 


In this section we shall show how arbitrarily small perturbations of the minimal model can 
induce chaotic dynamics of various types that are more substantial and complex than that 
described in Theorem 4.1. Our first perturbation involves the (geometric) local embedding 
of a Smale horseshoe in the minimal model in a neighborhood of the saddle point p*. 


5.1 Direct horseshoe chaos 


We begin by defining a map in a neighborhood of p* that yields the desired horseshoe. First, 
let 0 < (5o < 0.04 so that 

Be^Soip*) ■= {{x,y) E : {x-x^f+{y-y^f < 72So} C [a;*-6v^(5o, C (0,1)^ = int (/^) , 

where int(i?) denotes the interior of the subset E of the plane. Now, for each 0 < 5 < (5o,we 
introduce a (s, f)-coordinate system with origin at p*, the s-axis pointing to the right along 
the linear unstable manifold of p* and the t-axis pointing upward along the stable manifold 
of p*. Let Ps : {{s,t) : |s|, |t| < 65} —)■ be defined in terms of s, f-coordinates as 

Ps{s,t) := , (15) 

where 0^ is an odd function such that 

( -(3/2)s, 0 < s < 25 

(j)5{s) := < 3(s — 35), 25 < s < 45 , 

y (l/3)(s + 55), 45 < s < 65 

and 'ips is an odd function of t for each s given by 

ip5{s,t) := -St - 

where 

f 0, |s| < 25 

/i(s) ;= < sgn(s) (|s| — 5), 5 < s < 25 . 

[ sgn(s)5, |s| > 25 


As p* is a fixed point of both Ps and the minimal model F, for and given e > 0 we can 
choose 0 < 5 = 5(e) < 5o such that in terms of the euclidean norm, we have 


l|f (P) - Pi(p)ll < £ 
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for all p G -B 6 \/ 2 (So(p*)- Moreover, it is easy to see - as shown in Fig. [^- that the image of 
the square 

Qs := {{s,t) : —26 < s,t < 26} 


under the map (13), namely PsiQs), is a (double) horseshoe. Note that all of the above 
functions are continuous and piecewise linear, which means they have smooth approximations 
(obtained by smoothly rounding out the corners) that are arbitrarily C^-close to them. 
Consequently, we can and will assume that the perturbations chosen here to prove our next 
result and in the sequel are smooth. 


Theorem 6. For every e > 0 there exists a 0 < 6 = 6(e) < 0.04 such that the (smooth) 
perturbation of the minimal model map F defined as 


Fs(p) := (1 - p(r)) Ps(p) + p(r)F(p), 


where Ps is as in (15), r := ||p —p*|| and 


( 0, 0 < r < 46y/2 

p(r) := } (1/572) (r - 45^2) , 45^2 < r < 55^2 , 

\ 1, r > 5572 


(16) 


satisfies \\Fs(p) — F(p)\\ < e for all p E P and Fs (Qs) is a double horseshoe as shown in 
Fig. 0 


Proof. As noted above, the restriction 0 < 5 < 0.04 guarantees that Bg^siP*) contained 
in the interior of P and that by taking 5 sufficiently small, we can further insure that 
\\F(p) — P 3 (p)\\ < e on Hence, it follows from the dehnition of p and (16) that 

jjF 5 (p) — F(p)\\ < e on P. Finally, (16) implies that Fs (Qs) = Ps (Qs), which is the double 
horseshoe illustrated in Figj^ and this completes the proof. □ 


In light of Theorem the next result on the existence of horseshoe type chaos follows 
directly from the results of Birkhoff, Moser and Smale (cf. [I2l [THl ESI EEl ED], EE] 133]). 

Corollary 7. The perturbation Fs in Theorem^is chaotic on an invariant subset contained 
in the double horseshoe image described therein. In particular, Fs restricted to this invariant 
set is conjugate to the shift map on three symbols. 
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Figure 3: Embedded double horseshoe 

5.2 Snap-back repeller chaos 

There also are arbitrarily small perturbations of F exhibiting snap-back repeller chaos 

of the type described, for example, by Marotto [221123] . We start by showing how a small 
perturbation of F can turn p* into a source with at least four snap-back points and then 
show that it is possible to create inhnitely many snap-back points circling the hxed point. 
Once again we employ the s, f-coordinate system and 5o used in the preceding subsection 
to dehne the perturbation in a neighborhood of p*. In particular, we dehne Rs : {{s,t) : 
|s|, \t\ < 65} —)• as 

Rs{s,t) ■■= , (17) 

where (ps is an odd function of s dehned for s > 0 as 

(j)s{s) := -2s, 

and 'ips is an odd function of t dehned for f > 0 by 


'ipsit) := 


-2t, 0<t<6 

2{t-25), 5<t<36 , 

56 — t, 36 < t < 65 
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and graphed in Fig. 


As p* is a fixed point of both Rs and the minimal model F, for and given e > 0 we can 
choose 0 < 5 = (5(e) < (5o snch that in terms of the enclidean norm, we have 

||F(p)-/ii(p)|| <€ 

for all p G -B 6 F 2 ( 5 o(p*)- above, we are going to assnme with no loss of generality, that onr 
pertnrbations are actnally smooth, and this leads to onr next result, which is an analog of 
Theorem |6l 

Theorem 8. For every e > 0 there exists a 0 < 6 = (5(e) < 0.04 such that the {smooth) 
perturbation of the minimal model map F defined as 


rs{p) := (1 - ff(r)) Rs{p) + a{r)F{p), 


(18) 


where Rs is as in 


(17), r ■ = 


\\p - p^W and 


a{r) 


0 , 

(1/(572) (r - 5(572) , 

1 , 


0 < r < 5^72 
5(572 < r < 6(572 , 
r > 6^72 


satisfies \\Fs{p) — F(p)|| < e for all p G P and p* is a snap-back repeller for Fs having the 
chaotic dynamics described in IWf . 


Proof. The restriction 0 < (5 < 0.04 insures that Bq^sIp*) contained in the interior of P 
and that by taking 6 sufficiently small, we can further insure that ||F(p) — i? 5 (p)|| < e on 


P^VlsiP*)- Consequently, (T^) and the definition of the function a implies that ||F 5 (p) — F{p) || 
e on J^. It also follows from (18) that Fs{p) = Rs{p) for all p G B^^^^{p^), so p* is a hyper¬ 
bolic repeller for Fs having - as is clear from Fig. |^- four snap-back points at (s, t) = (0, ±25) 
and (s,t) = (0, ±55). The chaotic dynamics then follows from [22] and the proof is com¬ 
plete. □ 


< 


It should be noted that the chaos described in Theorem is essentially one-dimensional 
inasmuch as it is confined to the unstable manifold of p* with respect to the minimal model F. 
However, it is not difficult to see how the construction in the above proof can be modified to 
obtain higher dimensional snap-back repeller chaos. One need only consider a perturbation 
Rs given in polar form (with p* as the origin in the s, t- coordinate plane) as 

R 5 ■= 'f’si’i^) (cos 9, sin 9 ), 
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Figure 4: Coordinate function for snap-back repeller perturbation 


which has two full circles of snap-back points around p*. If in addition, we modify 'ips in the 
annulus 2(5 < r < 65 so that it is negative for certain ^-sectors, the location of the snap-back 
points near p* can be easily controlled. 


5.3 Chaos generated by embedding transverse homoclinic orbits 
and heteroclinic 2-cycles 


Chaos can also be generated by perturbing F so that it has transverse homoclinic points or 
transverse heteroclinic 2-cycle points, for then it exhibits chaotic subshift dynamics (c/. PQ 

[ini na [la EHi [33]). 

To begin, we prove a simple result showing how to create transverse intersections in a 
homoclinic or heteroclinic curve on a surface with an arbitrarily small perturbation. It 
should be noted that it is well known that such transverse intersections can be produced by 
arbitrarily small perturbations on general surfaces (c/. [ISIEHIES] ), but we shall find 
it useful for our simulations in the sequel to present a specialized result for the plane that 
is much easier to prove. In fact, the idea of the proof is quite transparent, involving just 
a carefully localized small sinusoidal perturbation normal to the homoclinic or heteroclinic 
curve (see Figj^, but the details are a bit involved. By a homoclinic curve or heteroclinic 
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curve joining a single, respectively, pair of distinct saddle points of a differentiable map of 
a differentiable surface, we mean an (open) curve contained in the stable (unstable) and 
(stable) manifolds of the single, respectively, pair of points that contains the saddle points in 
its closure. A heteroclinic 2-cycle is just pair of heteroclinic curves joining a pair of distinct 
saddle points. 

Lemma 9. Let f : S ^ S be a self-map of connected surface S in Mf. Suppose that 
p and q are identical or distinct saddle points of f joined, respectively, by a homoclinic or 
heteroclinic curve 7 that is contained both in hh“(p)nfh*(g) and a subset of S that is compact 
in If there is a point m_i G 7 such that uq := f{u-i), Ui := p{u-i), U 2 := P{u-i) and 
M3 := /"^(m-i) are contained in a connected open neighborhood U of the closed subarc k 0/7 
from «_! to U 3 and f is invertible on U, then there is an arbitrarily small perturbation g 
of f , equal to f except in an open subset V of U that contains a nontrivial closed subarc a of 
K from Uq to a point vq E j not containing ui, such that Wg{p) has a transverse intersection 
with Wg{q) in f{cr). 


Proof. First, we orient the curve 7 from p to g so that it follows the direction of positive 
iterates. Next, dehne the (oriented) arclength along 7 starting from uq to be s, so that 
s(mo) = 0 < s(mi) < 5 (^ 2 ) owing to the dehnition of 7 and the hypotheses. It follows from 
our assumptions that there exists a point u E n with 0 < s{u) < s(mi) such that the closed 
subarc x of 7 from uq to u satishes the following property: 


(PI) / ^(x), X) fix) and Pix) are pairwise disjoint closed subarcs of n. 


At each point u E j there is a unit normal vector that is unique if we specify an orientation, 
which we do by dehning the positive direction to be consistent with the right hand rule and 
the positive direction along the curve. We denote this vector by n{u), which now allows the 
dehnition of a positive and negative distance from 7 points x sufficiently close to 7 , which 
we denote by z/(x). We then combine this with a coordinate six) dehned to be s(0(a;)), 
where (fix) := u is the unique (for points sufficiently near 7 ) point on the curve such that 
the distance d(x, 7 ) from the point a; to 7 is equal to ||x — m||, where H-H is the Euclidean 
norm. Then basic results on normal bundles such as in [19] imply that there exists Aq > 0 
such that for every 0 < A < Aq, 

N\\= {x E S s(m_i) — A/2 < s(a;) < 5 (^ 3 ) + A/2; |z/(x)| < A} (19) 

is an open set contained in U and is a neighborhood of the closed subarc k, N\ does not 
contain p or g, and it is well dehned in the sense that every x E N\ is uniquely determined 
by its coordinates (s(a:), i/(a;)). We now have all the tools necessary to provide a simple 
description of small perturbations that possess the desired transverse intersections of 
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unstable and stable manifolds. In virtue of the differentiability of the map and the homoclinic 
or heteroclinic curve and the dehnition of Nx, there exists 0 < Ai < Aq such that if 0 < A < Ai 
and we dehne N\{x) ■= {x E S : 0 < s{x) < s{u); |z^(a:)| < A}, then the following obtains: 


(P2) / ^ [Nx{x)), Nx{x), f {Nx{x)) and p (-A^a(x)) are pairwise disjoint closed subsets of 
Nxo- 


Now for any 0 < e < Ai, we dehne the perturbation increment function : S' —)■ S' as 


A,(a;) := 


e (1 — A 


-1 




0 , 


X i Nxix) 


( 20 ) 


which is tantamount to saying that the s-coordinate function of A^ is zero and i/-coordinate 
function is zero except in Nxicr) where it is dehned to be the coefficient of the unit normal 


vector in (20). Note that this function also vanishes when s(a:) = 0, (1/2)s(m) or s(u); that 
is, all along the normals to the curve 7 at mq, h and a point tn G 7 with s{u) < s{w) = 


{l/2)s{u) < s(m). Moreover, the graph of ( 20 ) in the s, z/-plane has a transverse intersection 
with the s-axis at s = (l/2)s(-u). The desired perturbation of / is 


9e = f + A„ 


( 21 ) 


which is readily seen to have the desired properties. In particular, \\ge{x) — f{x)\\ < e for 
all X G S', Wpp) = Wpp) and Wpq) = Wpq) in a neighborhood of p and q, respectively, 
and the differentiability of both / and in a neighborhood of w together with (PI), (P 2 ) 
and the transversality property of A^ at w guarantee that there is a transverse intersection 
of Wpp) and Wpq) at f{w). Thus, the proof is complete. □ 



Figure 5: Transverse intersection perturbation 


As a hnal remark concerning the above lemma, it is rather easy to see how by smoothing 
corners and reducing the scale of the perturbation increment, if necessary, the perturbation 
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of the function can be chosen to be small. Moreover, it is a simple matter to extend 
the result to general surfaces using standard techniques from differential geometry and 
topology (c/. [IHl CH EOl I2HI1211 El])- 

Next we shall show how to embed C°-small dynamics in a neighborhood of the hxed point 
p* of the minimal model F that has a homoclinic orbit or a heteroclinic 2-cycle. Once done, 
we can apply Lemmato create chaotic dynamics in arbitrarily small perturbations of F. 
We shall use the time-one maps of the following Hamiltonian differential equations dehned 
for 5 > 0 as 


^ = p (^2 + ^2 ^ , 

V = - V^) , 


which has the Hamiltonian function 


Hs 


( 1 / 4 ) {e+i^y 


2 (e - 


( 22 ) 


(23) 


Here we have relabeled the s, t-coordinates used in subsections 5.1 and 5.2 as p-coordinates 
to avoid confusion with the time parameter t in (22). This equation yields a pair of homoclinic 
orbits corresponding to = 0 that are depicted in Fig. For the heteroclinic 2-cycle, we 
choose the system 


')=-/+ . (24) 


with Hamiltonian function 


lU = (1/4) [2 + ,f) - (x^lS'‘)] ■ 


(25) 


This system has a heteroclinic 2-cycle contained in = ^^/4, which is also shown in Fig. 
Now let the solution of (22) and (24) be denoted, respectively as 


and 


(^,h) = (^o,ho)) 

(26) 

(^,h) = it, (^o,ho)) • 

(27) 


The desired embeddings are obtained from the time-1 maps of the above; namely, we 
dehne 


^5ii,'n) ■= Vi (1, (^,h)) 

(28) 

■='iPiihii,v)) ■ 

(29) 
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and 







Clearly, these maps have the desired homoclinic and heteroclinic orbits, respectively. Now, 
we select 5o > 0 so small that Bqsq{p^) is contained in the interior of P and then 0 < 5 = 
^(e) < ^0 for a given e > 0 snch that 



\\F{z)-^s{z)\\<e 

(30) 

and 

\\F{z)-4>s{z)\\<e 

(31) 

for all z e Basoip*)- 

Whence, we can dehne 



Fs{z) := (1 - x(r)) 4 ) 5 ( 2 ;) -h x{r)F{z) 

(32) 

and 

Fs{z) := (1 - x(r)) 4 ) 5 ( 2 ;) -F x(r)F( 2 ;), 

(33) 

where r ■= \\z — p*| 

and 



x(r) := < (1/(5) (r — 3S ), 36 < r < 46 


(34) 


1, r > 46 




Figure 6: Embedding homoclinic and heteroclinic orbits 


It follows from (24)-(33) that 
and that Fs has a 


F(z) - Fi(z) 


< e and 


F(z) - h(z) 


< e for all z e B 


lomoclinic orbit comprised of the stable and unstable manifolds at p* 
and Fs has a heteroclinic 2-cycle centered at p* as shown in Fig. Whence, we may, 
after choosing 6 smaller if necessary, use Lemma ^ to further perturb Fs and Fs to create a 
transverse intersection point in the homoclinic orbit and, respectively, to create a transverse 
intersection point in one or both components of the heteroclinic 2-cycle, while still remaining 
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e-close to the original maps in the (7° metric. Then A-Lemma based arguments of the type 
developed in such sources as P [ini [281 [33] and the smoothing discussed above lead directly 
to the following result. 

Theorem 10. For every e > 0 there exist smooth e-close perturbations of the minimal 
model map F exhibiting transverse homoclinic point or transverse heteroclinic 2-cycle induced 
chaos. 


5.4 Multihorseshoe strange attractor perturbations 

The next element of our description and analysis of more pervasive chaotic perturbations 
of the map F dehned by 0-0 involves embedding a symmetric pair of attracting horseshoes 
along the lines introduced in [T^. In particular, we shall show how to perturb F to 
produce multihorseshoe (more precisely, double-horseshoe) chaos (c/.fTH], as shown in Fig. 

m 


To begin the construction of the embedding, we note hrst that the minimal map F 
is actually dehned and smooth at all points in the x, y-plane above the curve dehned by 
X y — xy = t). Next, for convenience, we denote the period-2 points (1, 0) and (0,1) by p 
and g, respectively. Let e > 0 be given. As both p and q are hxed points of F^ and of 
where G is the rehection in the line y = x, we may choose 0 < 5 < 1/4 such that 


||F(;.)-G(;.)||<e/2 

whenever ^ G U Dehne 

f ^ ^ ^10725(P) ^ ^10725(?) 

F,{z) := < (1 - uj{rp)) G{z) + uj{rp)F{z), z G B^o^^ip) 

( {^-^irq))G{z)+u:{rq)F{z), z e B^f^^^{q) 

where rp := \\z — p\\, rq := \\z — g|| and 

0, 0 < r < 85^ 

uj{r) := <{ (1/5^2) (r - 85^2) , 85^2 < r < 95^2 . 

1, r > %y/2 


It is clear from the dehnition and (35) that 


F{z)-F,{z) <6/2. 


(35) 


(36) 


(37) 


(38) 
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It now remains to make a final adjustment of the map that has the attracting horse¬ 
shoes. To this end, we dehne the following perturbation of the identity in n P 

using the original x, ^/-coordinates 




where 


and 


ip{x) 


1 , 

2{x -F 5) - 1, 
-2{x + 75) + 3, 


a:<l — 75orl — (5<x<l 
1 — 4:5 <x<l — 5 
l — 75<x<l — 45 




I//5, 

{y/7>) -h6(l -5(5-x), 
{y/h) -F 65 

(///5) + 3 (a: + 95 - 1) 


a;<l — 95orl — 55<a;<l 
1 — 65<a;<l — 55 
1 — 75<a;<l — 65 
1 - 95 < X < 1 - 75 


(39) 


(40) 


Observe that 0 has a sink at p and a saddle point at (1 — 25, 0) with a horizontal unstable 
and vertical stable manifold. Moreover, 0 maps the rectangle [1 — 45,1 — 5] x [—5, 95] onto 
a piecewise smooth (and of course smoothable) attracting horseshoe as dehned in [15]. This 
function can be reflected in the line y = x to obtain the symmetric map 0 := G o Q o G 
mapping the rectangle [—5, 95] x [1 — 45,1 — 5] onto the reflection of the horseshoe image of 
Q in y = X. 


Now it follows directly from the above dehnitions that, taking 5 smaller if necessary, we 
can insure that 

||G(^)-0oG(z)|| <e/2 (41) 

when z G fl P and 


G{z)-QoG{z) 


< e/2 


(42) 


for all z G i?iov^ 5 (p) 0 P- Therefore, also taking into account 
dehned as 


(38), the modihcation of T). 


[ ^ ^ ^ ^10V2<5(p) U ^10^25(?) 

^e(^) := < {l-uj{rp))QoG{z)Fuj{rp)F{z), z e 

I {l-i^{rq))QoG{z)Fuj{rq)F{z), z ^ 

satishes 

||F(x)-F.(x)||<e. 

Moreover, we note that both p and q are attracting hxed points (sinks) of F/, just as they 
are of F^, while (1 — 25, 0) and (0,1 — 25) comprise a 2-cycle of F^ (consisting of saddle points 
of F/) but not F. Finally, as constructed, each of the symmetric horseshoes shown in Fig. 
E is an attracting horseshoe of F, so that it follows from the main multihorseshoe theorem 
in [T5| that we have now proved the following result. 
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Theorem 11. For every positive e there is a smooth e - close perturbation of the min¬ 
imal map F having a strange chaotic double-horseshoe attractor with attracting horseshoes 
in neighborhoods of the points (1,0) and (0,1). 


It is interesting to note that the strange attractor for the small pertnrbation of F in 
Theorem resembles a discrete analog of the “donble scroll” attractor for Chna’s circnit 
{cf [6l[7]), also seen in the dynamics of RSFF realization simnlations snch as in |1] and 
conld, with some minor modihcation, produce discrete analogs of the attractors found in the 
various physical circuit (ODE) models such as in [121 HZl EH ESI EZj- 



Figure 7: Embedded double-horseshoe chaotic strange attractor 


5.5 Neimark—Sacker bifurcation perturbations 


Our hnal result shall show that Neimark-Sacker bifurcations can occur in arbitrarily small 
perturbations of the minimal model map F. Once again, we use the r^-coordinates 
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employed in preceding subsections to define for 5 > 0 the following parameter-dependent 
map in polar coordinates in the - plane with origin at p* 


Ns (^, p; /i, a) := —S tanh(/ir/5) (cos(6' -f- a), sin(6* -|- a)), 


(43) 


where 1/2 < /i < 3/2, a is nonnegative and r := ||(^, p) — P*||- It is easy to verify that the 


origin is a global attractor of (43) for 1/2 < /i < 1 and a local repeller for 1 < p < 3/2, 
so /i = 1 is a bifurcation value. Moreover, as fi increases across 1 with a corresponding 
transition of the origin from a sink to a source, a stable invariant circle of radius r = r(/r), 
where r(/r) is the unique positive solution of 


6 tanh(pr/5) = r. 


(44) 


The parameter a just represents the rotation of the map (43), so what we have is a Neimark- 
Sacker bifurcation at /r = 1. 


Now just as in the preceding subsections, we can choose > 0 so small that Bqsq{p*) is 
contained in the interior of N and then 0 < 5 = ^(e) < for a given e > 0 such that 


\\F{z) - Ns{z)\\ < e 

for all fi G (1/2, 3/2) and a > 0 whenever G BqSq{p*)- Therefore, by defining 

Fs{z) := (1 - x(r)) Ns{z) + x{r)F{z), 


(45) 


where x is defined just as in (34), we obtain a map that is e-close to F in the norm and 


has the desired bifurcation properties. In short, we have now proved the following result. 

Theorem 12. For every positive e there is a smooth e - close perturbation F of the 
minimal map F having a Niemark-Sacker bifurcation at p*. 


A direct detailed construction was used for the proof of Theorem but the same result 
can be proved, using the main theorem in [5] , for any perturbation in which p* changes from 
an attractor to a local repeller as a parameter is varied. It should also be noted that one 
can, by a straightforward modification of the above procedure, construct arbitrarily small 
(7° perturbations of F exhibiting doubling Neimark-Sacker (Hopf) cascades like those in the 
ad hoc RSFF model analyzed in [2]. 


6 Simulations, Computations and Comparisons 


Our purpose in this section is to show with just a few examples that our discrete dynamical 
model - when properly perturbed - shares many properties with actual physical realizations 
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(and their associated mathematical models) of the R-S flip-flop and related circuits such as 
in [HI Uni [TTl [211 IMl |26l [271 jHH [3l]. For example, we have already demonstrated in Section 
4 that our model can be perturbed so that it exhibits the chaos found by simulation in a 
one-dimensional map associated with the dynamics of the realization of the flip-flop circuit 
- which is not an R-S flip-flop circuit - investigated in [27], and it is this kind of chaos we 
consider in the next subsection. 


6.1 Perturbed minimal model with one-dimensional chaos 


We derive the perturbed map by plugging in 


(13) to (|^ where x —>■ (p^ix) and y —>■ (psiy), 


^{x,y) 

r]{x,y) 


4>5{y){l - (pajx)) 

1 - {1 - (p^{x)){l - (ps{y)) 
(pa{x){l - (psjy)) 

1 - (1 - (pa{x)){l - (psiy))' 


(46) 


The iterates of the perturbed map are shown in Fig. 8b 


Iteration Number = 50 



Xn 


(a) Illustration of Eq. 


13 


(b) Plot of the iterates of the perturbed map. 


Corollary 13. The perturbation in the preceding theorem can be chosen so that it is a C° 
function having the same gualitative dynamics as the function constructed above. 


6.2 Evidence of two-dimensional chaos 


As we have shown above, almost any type of chaotic dynamics - including double scroll 
chaos - can be obtained from the ideal map by inserting specihc localized perturbations. We 
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shall now show that the introduction of a fairly general type of small (7° perturbation is 
apt to produce chaotic dynamics. In order to produce two-dimensional chaos lets use the 
perturbation, 


k 

= 3.7[ao — .02 — ajCosjTTx] 


^0,cr (^) 
(^) 


f^x ^x I 2 ^ 

8 ^ 4 

^{(1 - A^)sinj(^ - e)7r + [(-1)^ 

ZJTl 2 J7l 


1 ]} 


to get the map, 


1 - (1 - ipm,5{y))i^ - 
(p^,,(x)(l - iPm,s{y)) 
1 - (1 - ipm,5{y))i^ - 


(47) 


The iterates of this map, with their characteristic splattering indicative of chaos, are 
shown in Fig. 


Iteration Number = 1000 



Figure 8: Iterates of a perturbed map seeming to exhibit two-dimensional chaos. 
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6.3 Ring oscillator example 


Another interesting circuit, which we intend to analyze in detail in a forthcoming paper, is 
a modihed ring oscillator. This ring oscillator comprises three NOR gates with feedbacks, 
in a very similar fashion to that of the RS flip-flop circuit. One may even choose to think of 
this as a three-dimensional RS flip-flop circuit. 



Figure 9: A schematic of a ring oscillator designed out of three NOR gates. 

We note that if the inputs are set to zero, the system is precisely a ring oscillator because 
the NOR gates now act as inverters. Applying our algorithm for hnding a discrete dynamical 
system model for logical circuits yields the following “ideal” (unperturbed) model. 


r] = 

c = 


{x - l){y{z - 1) + 1) 

(z-l)(x(?/-1)-hi) 


(48) 


The only valid fixed point is = ((3 — ■\/5)/2, (3 — a/5)/2, (3 — a/5)/2). Since 

the fixed points are roots of a quartic equation, one may wonder what happens to the other 
roots. We hnd that the other three roots are out of our domain, which means that they 
are of no practical consequence. For the ideal model, regardless of the initial conditions, the 


orbits decay to the hxed point in a spiral manner shown in Fig. 10 


In order to produce more interesting dynamics we perturb the ideal model slightly. Using 


a perturbation similar to (47), we produce the chaotic dynamics shown in Fig. 11 













»r = 20 


9r = S3 



Figure 10: Plots of the orbits of Eq. 
(.1, .1, .9) respectively. 


48 with initial conditions, (.99, .1,.9), (.1,.9, .2), and 


Iteration N urn ber= 100000 



Figure 11: Plots of the orbits of the perturbed model with initial conditions (.99, .1, .9). 
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7 Concluding Remarks 


We have introduced and analyzed a rather simple discrete (ideal) dynamical model - grounded 
on hrst principles - for the RSFF circuit, which is based upon the iterates of a planar map. 
Moreover, we have proved that this model can be modihed - by arbitrarily small pertur¬ 
bations - to produce just about any dynamical property observed in physical realizations of 
RSFF and related flip-flop circuits. We have also shown that rather general small pertur¬ 
bations are apt to change the structurally stable ideal model to a two-dimensional discrete 
dynamical system with chaotic dynamics and related artifacts such as strange chaotic at¬ 
tractors. 

Naturally, we are planning to extend this discrete dynamical systems approach to a much 
wider class of logical circuits and their perturbations, and verify the effectiveness of this ap¬ 
proach by comparing our dynamic predictions with those of more standard ODE approaches 
as well as experimental data extracted from actual physical circuit measurements. We also 
plan to address a number of related questions such as showing that logical circuit realizations 
represent perturbations of our ideal models in some sense, and that the perturbations can 
actually be characterized and quantihed. Such an investigation should provide insight into 
the relationship that we believe exists between our discrete dynamical systems approach and 
the underlying discrete dynamics of reconhgurable chaotic logic gates [Ill- 
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